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We show that the first year results of the Wilkinson Microwave Anisotropy Probe (WMAP) 
constrain very efficiently the energy density in relativistic particles in the universe. We derive new 
bounds on additional relativistic degrees of freedom expressed in terms of an excess in the effective 
number of light neutrinos AN e s. Within the flat ACDM scenario, the allowed range is AN c h < 6 
(95% confidence level) using WMAP data only, or -2.6 < AN cfi < 4 with the prior H = 72 ± 8 
km s _1 Mpc -1 . When other cosmic microwave background and large scale structure experiments 
are taken into account, the window shrinks to —1.6 < AN c g < 3.8. These results are in perfect 
agreement with the bounds from primordial nucleosynthesis. Non-minimal cosmological models with 
extra relativistic degrees of freedom are now severely restricted. 

PACS numbers: 14.60.Pq, 98.70.Vc, 98.80.Es 



> 

co 
m 

(N 
O 
CO 

O 

Or 

6 ■ 
a: 

C/3 . 



I. INTRODUCTION 

What is the matter budget of the universe? This fas- 
cinating and central question in cosmology is currently 
being answered with increasing precision, thanks to out- 
standing measurements of the cosmic microwave back- 
ground (CMB) anisotropies, correlated with the study of 
large scale structure (LSS), of the primordial abundances 
of light elements, and many other observables in the far 
universe. 

One of the most intriguing issues is to determine the 
contributions of matter and dark energy to the total en- 
ergy density (fl m and fl\ in units of the critical density). 
Inflation predicts Q m + Q\ = 1 but, so far, there is no 
theoretical prediction on the value of each of these pa- 
rameters. This explains why in most cosmological pa- 
rameter analyses they receive much more attention than 
the radiation density f2 r , which is often assumed to be 
well-known. However, the reference value of fl r relies 
on a strong theoretical prejudice: apart from the CMB 
photons, the dominant relativistic backgound would con- 
sist of the three families of neutrinos, whose temperature 
would be fixed with respect to the CMB temperature by 
the standard picture of neutrino decoupling prior to Big 
Bang Nucleosynthesis (BBN). 

Measuring the energy density of radiation today is not 
easy, because it is known to be three orders of magnitude 
below the critical density. The main constraints come ei- 
ther from the very early universe, when radiation was the 
dominant source of energy, or from the observation of cos- 
mological perturbations, which imprint the time of equal- 
ity between matter and radiation. In particular, through 
BBN models, the primordial abundances of light elements 
can be related to f2 r , evaluated at the time when the 
mean energy in the universe was of order 1 MeV 4 - while 
the power spectrum of photon anisotropies and of matter 
density carry a clear signature of the time at which fl r 
became comparable to Q m , at energies of order (0.1 eV) 4 . 



Until now, the BBN constraint on fl r was more stringent 
than the one from cosmological perturbations, by approx- 
imately one order of magnitude. This leaves the door 
wide open for various plausible assumptions concerning 
the radiation content of the universe, which should not 
be necessarily the same during BBN and at the time of 
matter/radiation equality. For instance, a population of 
non-relativistic particles may decay into relativistic ones, 
enhancing the radiation energy density. Moreover, the 
standard BBN scenario itself - which is the simplest way 
to explain the formation of light elements, but not the 
only one - needs to be tested. In this respect, the best 
would be to have some independent measurements of the 
two free parameters of BBN: the baryon density f^ft 2 
and the radiation density Q, r h 2 (where the reduced Hub- 
ble constant is h = H /(100 km s _1 Mpc -1 )). The later 
controls the expansion rate in the early universe. Some 
precise bounds on fl^h 2 were already obtained from re- 
cent CMB experiments, while the determination of Q, r h 2 
was still quite loose, compared to BBN predictions. 

The goal of this paper is to update this analysis and 
to show that the outstanding data from the first year 
sky survey of the Wilkinson Microwave Anisotropy Probe 
(WMAP) [liQ gives better constraints, showing increas- 
ing evidence in favor of standard BBN, and leaving very 
small room for extra relativistic degrees of freedom be- 
yond the three neutrino flavors. 



II. THE EFFECTIVE NUMBER OF 
RELATIVISTIC NEUTRINOS 



The energy density stored in relativistic species, p r , 
is customarily given in terms of the so-called effective 
number of relativistic neutrino species N e g (see [3j for a 
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review and references), through the relation 
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where p 7 is the energy density of photons, whose value 
today is known from the measurement of the CMB tem- 
perature. Eq. can be also written as 
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where p% denotes the energy density of a single species of 
massless neutrino with an equilibrium Fcrmi-Dirac distri- 
bution with zero chemical potential, and p^ is the pho- 
ton energy density in the approximation of instantaneous 
neutrino decoupling. The normalization of iV e fj is such 
that it gives iV e g = 3 in the standard case of three fla- 
vors of massless neutrinos, again in the limit of instanta- 
neous decoupling. In principle N B s includes, in addition 
to the standard neutrinos, a potential contribution p x 
from other relativistic relics such as majorons or sterile 
neutrinos. 

It turns out that even in the standard case of three neu- 
trino flavors the effective number of relativistic neutrino 
species is not exactly 3. The decoupling of neutrinos from 
the rest of the primordial plasma occurs at a temperature 
of 2 — 3 MeV, not far from temperatures of order the elec- 
tron mass at which electron-positron annihilations trans- 
fer their entropy into photons, causing the well-known 
difference between the temperatures of relic photons and 
relic neutrinos, T/T v = (H/4) 1 ^ 3 (see e.g. 0). Accurate 
calculations [1H 0] nave shown that neutrinos are still 
slightly interacting with e ± , thus sharing a small part of 
the entropy release. This causes a momentum dependent 
distortion in the neutrino spectra from the equilibrium 
Fermi-Dirac behavior and a slightly smaller T /T v ratio. 
Both effects lead to a value of N e g = 3.034. A further, 
though smaller, effect on T/T u is induced by finite tem- 
perature Quantum Electrodynamics (QED) corrections 
to the electromagnetic plasma [H, Q • A recent combined 
study of the incomplete neutrino decoupling and QED 
corrections concluded that the total effect corresponds 
to N cS = 3.0395 ~ 3.04 Therefore we define the 

extra energy density in radiation form as 
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The standard value of N e s corresponds to the case of 
massless or very light neutrinos, i.e. those with masses 
much smaller than 1 eV. More massive neutrinos affect 
the late avolution of the universe in a way that can not 
be parametrized with a AN c g . However, the recent evi- 
dences of flavor neutrino oscillations in atmospheric and 
solar neutrinos, in particular after the recent KamLAND 
data show that the neutrino masses are not large enough, 
except in the case when the three mass eigenstates are 
degenerate (see e.g. for a recent review). We do not 



consider such a case in the present paper, but assume 
that the neutrino mass scheme is hierarchical, with the 
largest mass of order m„ ~ yj Am^ tm ~ 0.05 eV. 

The value of AN c s is constrained at the BBN epoch 
from the comparison of theoretical predictions and exper- 
imental data on the primordial abundances of light ele- 
ments. Typically, the BBN bounds are of order AN e g < 
0.4-1 El Ellll E El. Independent bounds on the 
radiation content of the universe can be extracted from 
the analysis of the power spectrum of CMB anisotropies. 
An enhanced contribution of relativistic particles delays 
the epoch of matter-radiation equality, which in turn in- 
creases the early integrated Sachs- Wolfe effect. Basically 
this leads to more power around the scale of the first 
CMB pea k. Previous analyses found weak bounds on 
AAeff [H HI [H HI, HU , that can be significantly im- 
proved by adding priors on the age of the universe or by 
including supernovae and LSS data (2^ • One of the most 
recent bounds on N e s, from a combination of CMB and 
PSCz data [U, is N cS = 6±f. 5 (95% CL). Thus these 
bounds were not as restrictive as those from BBN. How- 
ever, the precise measurements of WMAP (and those of 
PLANCK in the near future) are going to significantly 
improve the CMB constraint on AN c g, as shown in var- 
ious forecast analyses (see for instance |23|, |24, 25] ) and 
in the calculations of this paper (see section ITV|) . 

Many extensions of the Standard Model of particle 
physics predict additional relativistic degrees of freedom 
that will contribute to AN c g. There exist models with 
4 neutrinos which include an additional sterile neutrino 
in order to explain the third experimental indication of 
neutrino oscillations (the LSND results) . It was shown 
in many studies (see for instance [HEll) that a ^ f° ur 
neutrino models, both of 2+2 and 3+1 type, lead to a full 
thermalization of the sterile neutrino flavor before BBN, 
and thus to AN c g ~ 1, a value disfavored in the standard 
minimal model of BBN. Moreover, in these models there 
exists at least one neutrino state with mass of order 1 
eV. 

It is also possible that the relativistic degrees of free- 
dom at the BBN and CMB epochs differ, for instance 
because of particle decays which increase the photon tem- 
perature relative to the neutrino one [2*sj. In some situ- 
ations AA^ff can be effectively negative at BBN, such as 
the case of a distortion in the v P or v e spectra [2! |3(J , 
or a very low reheating scenario |31| . 

A non-standard case that has been considered many 
times in the past is the existence of relic neutrino asym- 
metries, namely when the number of neutrinos and an- 
tineutrinos of the same flavor is significantly different. 
These so-called degenerate neutrinos are described by a 
dimensionless chemical potential £ a = fj, Ua /T, and it has 
been shown that the neutrino energy density always in- 
creases for any value £ Q ^ 
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Interestingly, some combinations of pairs (£e,£^,i") could 
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still produce the primordial abundances of light elements 
for a larger baryon asymmetry, in the so-called degen- 
erate BBN scenario At the same ti me, the weaker 
CMB bounds on are flavor blind |24l l33|. However, 
it was recently shown that for neutrino oscillation pa- 
rameters in the regions favored by atmospheric and solar 
neutrino data flavor equilibrium between all active neu- 
trin o sp ecies is established well before the BBN epoch 
[HHEHil. Thus the stringent BBN bounds on £ e ap- 
ply to all flavors, so that the contribution of a potential 
relic neutrino asymmetry to AiV e g is limited to very low 
values. 



III. METHOD 

Using the cmbfast code [37j, we computed the cos- 
mological perturbations (temperature and polarization 

anisotropies c\ T ' E '' 1 , matter power spectrum P{k)) for 
a grid of models with the following parameters: baryon 
density Ub = Qbh 2 , cold dark matter density cu c dm — 
Qcdmh 2 , hubble parameter h, scalar tilt n s , optical depth 
to reionization r, global normalization which is not dis- 
cretized - and of course an additional contribution from 
relativistic particles AN e g. We include corrections to 
the CMB spectra from gravitational lensing [3^ > as com- 
puted by cmbfast. Apart from AiV c ff , our set of param- 
eters is the simplest one used by the WMAP team in 
their parameter analysis Q, and accounts very well for 
the first year WMAP data. We restrict ourselves to a 
flat universe: since the curvature is known to be small 
from the position of the first CMB peak, we adopt the 
theoretical prejudice that the universe is exactly flat, as 
predicted by inflation, rather than almost flat. There- 
fore, fl\ is equal to 1 — (ajt> + ^cdm)/h 2 . Allowing for 
a small curvature could alter our results by a few per- 
cent. We also neglect the possible contribution of grav- 
itational waves and a possible scale-dependent tilt. A 
running tilt, in favor of which the WMAP collaboration 
finds some marginal evidence, would not change our pre- 
dictions based on WMAP alone, because the later does 
not constrain the primordial spectrum on a wide enough 
range of scales. However, it could slightly alter our re- 
sults based on CMB and LSS data. 

Our grid covers the following ranges: 0.019 < Mb < 
0.028, 0.065 < uj cdm < 0.27, 0.5 < h < 0.9, 0.8 < n s < 
1.28, < r < 0.5, -3 < AiV cff < 5. We analyze it using 
an interpolation and minimization routine developed at 
LAPTH. Our code performs a multi-dimensional interpo- 
lation for each value of Ci or P(k) in order to obtain the 
spectrum at any arbitrary point, and then, computes the 
likelihood of the model. For WMAP, the likelihood is cal- 
culated using the software kindly provided at the NASA 
web site |39j, and explained in |40j. We will also define 
a combined likelihood including the pre-WMAP CMB 
data compilation by Wang et al. |4l| (which is still use- 
ful for constraining high multipoles), and the LSS data 
derived by Percival et al. |42] (32 points on wavenumbers 
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FIG. 1: The A7V cff likelihood for WMAP + weak h prior 
(green), WMAP + strong h prior (red), and the same plus 
other CMB experiments and the 2dF redshift survey (blue). 
The horizontal lines show the 68% (resp. 95%) confidence lev- 
els. The step of AiV e ff in our grid of models is 0.5. 

k < 0.15 Wpc^ 1 ) from the 2dF redshift survey Hjj. For 
these two data sets, we use window functions and corre- 
lation matrices available from |4*3.l45j . We constrain each 
free parameter using a Bayesian approach: the 68% (resp. 
95%) confidence limits are defined as the values for which 
the marginalized likelihood drops by exp[— (x 2 , — l)/2] 
(resp. exp[— (xo — 4)/2]), where Xo is the best chi-square 
value in the whole parameter space. 



IV. RESULTS AND DISCUSSION 

We start the analysis using only the first year WMAP 
temperature and polarization data, plus a weak prior 
0.5 < h < 0.9, which is implicit from the limita- 
tion of the grid. We checked that for AN c g = 0, 
we find the same bounds as the WMAP collaboration, 
with a minimal effective chi square XcS ~ 1431.5 for 
1342 effective degrees of freedom (d.o.f.). The best 
fit over our whole grid has (AN e s,u>b,h,ui c d m ,T,n s ) = 
(-0.9,0.024,0.68,0.10,0.16,0.97), and a x 2 value which 
is not significantly lower (XcS — 1431.1 for 1341 d.o.f.), 
showing that a non-zero A.N e g is not required in or- 
der to improve the goodness of fit. At the 95% confi- 
dence level (CL), we can derive only an upper bound 
AiVeff < 6, which is impressively smaller than the pre- 
vious bound from CMB only, AiV cff < 14 HIHH. At 
the l-o" level, the extra relativistic energy density is lim- 
ited to the range —2.2 < AiV c ff < 1.2, corresponding 
to a dispersion of 3.4. This is in nice agreement with 
the prediction from ^ that with the full WMAP data, 
one would reach a dispersion of 3.17. The analysis re- 
veals that the indetermination of this parameter is caused 
mainly by a degeneracy with h, as shown in previous 
works. When AiV e ff runs from -3 to 6, the best-fit value 
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FIG. 2: The two-dimensional confidence limits on (AN e g, 
h), based on CMB and LSS data, at the l-cr (dark blue) and 
2-cr (light blue) levels. The superimposed yellow stripe shows 
the HST Key Project measurement of h (l-cr level). 



TABLE I: The best fit values and 2-cr (95% CL) limits on 
AiVcfi for the different data sets. 



data set 


AN cS 


WMAP + weak h prior 


-0.9t!;5 


WMAP + strong h prior 


-o.Tifi 


WMAP + other CMB + LSS + strong h prior 





of ft, goes from 0.55 to 0.90 (while oj c d m decreases, main- 
taining an approximately constant value of the matter 
fraction Q m = 0.31 ± 0.03). After AN cS = 6, the weak 
prior h < 0.9 is saturated, and the probability drops 
abruptly. 

In order to remove the degeneracy, we add to the def- 
inition of the effective \ 2 a gaussian prior on h derived 
from the HST Key Project gfj, h = 0.72 ± 0.08 (at the 
l-cr level). This prior is sufficient to reduce significantly 
the 95% allowed window: -2.6 < AN eS < 4.0. The no 
neutrino case with AN c s = — 3 is compatible with the 
data only beyond 99% CL. The best-fit model has still 
a slightly negative AN c g = —0.7, but this is not statis- 
tically significant. The best xls does not change very 
much (1431.1 for 1342 d.o.f.) because WMAP alone is in 
remarkable agreement with the Key Project value [2|. 

Finally, when we include the pre- WMAP data compila- 
tion and the 2dF redshift survey (i.e., information on the 
third and fourth CMB peak, and on the scale of the turn- 



over in the matter power spectrum), the window tends 
to shift a little bit towards positive values of AA^ff . The 
best-fit model has AN cS = 0.5 and xls = 1492 for 1396 
d.o.f. The 95% allowed window is -1.6 < AAr cff < 3.8. 
Using all data, we find that the no neutrino case is ex- 
cluded at 99.9% CL, which constitutes a clear indication 
of the presence of relic background neutrinos, already 
shown by pre- WMAP data Hcj. 

Note that throughout the analysis, we did not include 
any constraint from supernovae data. This is because all 
good-fitting models have naturally (fi m , ^a) very close 
to (0.3, 0.7). So, adding a supernovae prior would be 
completely irrelevant - unlike the h prior, which plays 
a crucial role in removing the degeneracy with AA c ff. 
In order to emphasize this last point, we show in fig- 
ure [21 the two-dimensional l-cr and 2-cr confidence limits 
in (AJVcs, h) parameter space. The contours are based 
only on the CMB and LSS data (the HST Key Project 
result is just superimposed as a yellow stripe). The de- 
generacy between these two parameters clearly appears, 
showing that any improvement in the direct determina- 
tion of h will be crucial for closing the AiV c g allowed 
window. For instance, large deviations from the standard 
value AA c ff = would be excluded if h were measured to 
be very close to 0.70, while h > 0.75 would bring strong 
evidence for extra relativistic species. 

To summarize, we show in table [I] the best fit values 
and 2-cr (95% CL) limits on AA^g for the three different 
data sets. Our results show that the first year results 
of WMAP significantly improve the bounds on an 
additional contribution to the radiation density of the 
universe. In the near future, the updated data from 
WMAP and the new results from the Sloan Digital 
Sky Survey 01 should allow for an even better deter- 
mination, providing a high-precision test of primordial 
nucleosynthesis . 

Note added: After the submission of this work, our 
results were confirmed in refs. (with a generalization 
of the bounds to a non-flat Universe) and [43 • 
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